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A search for the pair production of scalar top quarks, t, has been performed in 360 pb -1 of data 
from pp collisions at a center-of-mass energy of 1.96 TeV, collected by the DO detector at the Fermilab 
Tevatron collider. The t decay mode considered is t — > cxi, where x? is the lightest supersymmetric 
particle. The topology analyzed therefore consists of a pair of acoplanar heavy-flavor jets with 
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missing transverse energy. The data and standard model expectation are in agreement, and a 95% 
C.L. exclusion domain in the (m^m^o) plane has been determined, extending the domain excluded 
by previous experiments. 

PACS numbers: 14.80.Ly, 12.60.Jv 



Supersymmetric (SUSY) models [l| predict the exis- 
tence of new particles, carrying the same quantum num- 
bers as their standard model (SM) partners, but differing 
by half a unit of spin. For instance, there are two scalar- 
quark fields associated with the left- and right-handed de- 
grees of freedom of each ordinary quark. The mass eigen- 
states result from the diagonalization of a mass matrix, 
with elements determined by the specific SUSY-breaking 
pattern. A light SUSY partner of the top quark, or stop, 
is a generic prediction of models in which the scalar quark 
masses are equal at the grand unification scale. A first 
reason is that, due to the impact of the large top quark 
Yukawa coupling in the renormalization group equations, 
the diagonal elements of the mass matrix are driven to 
values smaller than those for the other scalar quarks at 
the electroweak scale [2J. A second reason is that the 
off-diagonal terms are proportional to the relevant quark 
mass, and hence are much larger in the case of the top 
quark. The mass eigenstates are therefore broadly split, 
with the mass of the lighter stop t thus driven to an even 
lower value [3j . Finally, a light stop is a necessary ingre- 
dient in the context of electroweak baryogenesis [J] . 

In models with i?-parity conservation [5j] , the lightest 
SUSY particle (LSP) is stable, and cosmological con- 
straints imply that it should be neutral and colorless Q . 
In a large class of SUSY models, the lightest of the neu- 
tralinos — the mass eigenstates resulting from the mix- 
ing of the SUSY partners of the neutral gauge and Higgs 
bosons — is the LSP, which furthermore appears as a 
viable dark matter candidate. In the following, it will be 
assumed that i?-parity is conserved and that the LSP is 
the lightest neutralino Xv 

The dominant stop decay modes are expected to be 
i ~ * tXi and t — > bxi , where the chargino is the 
lighter of the two mass eigenstates resulting from the 
mixing of the SUSY partners of the charged gauge and 
Higgs bosons. However, in the t mass range of interest in 
this Letter, the t — > tx® decay mode is kinematically for- 
bidden. In the following, the region of SUSY parameter 
space with < nib + m ~+ and < Mw + nib + m^o 
is considered, and it is assumed that t — > ex? , a flavor- 
changing loop decay 0], is the only relevant decay mode, 
i.e., that the tree- level four-body decays [§| t — ► bff'xi 
can be neglected. 

In pp collisions, stop pair production proceeds via qq 
annihilation and gluon-gluon fusion. The cross section 
has very little dependence on SUSY parameters other 
than the stop mass. At the center-of-mass energy of 
1.96 TeV available in Run II of the Fcrmilab Tevatron col- 



lider, it ranges from 15 to 2.25 pb for stop masses from 
100 to 140 GeV, as calculated at next-to-leading order 
(NLO) with PROSPINO 0, for equal renormalization and 
factorization scales fj, r f = mi and using the CTEQ6.1M 
parton distribution functions (PDFs) Jl(| . The final state 
topology resulting from the t — > cxi decay is a pair of 
acoplanar jets, with large missing transverse energy JpT 
carried away by the two weakly interacting LSPs. Previ- 
ous searches in this topology performed at LEP excluded 
stop masses smaller than k, 100 GeV, essentially inde- 
pendent of the stop-Xj 1 mass differenc e 111! . Searches in 
data from the Run I of the Tevatron [l2l [l3| extended 
the domain excluded at LEP to larger stop masses, but 
for Xi masses not exceeding « 50 GeV. The largest stop 
mass excluded was 122 GeV, for m^o — 45GeV[l3j]. In 
this Letter, we report on a similar search, performed in 
data collected using the DO detector during Run II of the 
Tevatron. 



The acoplanar jet topology may arise from new physics 
processes other than stop pair production. Recently, the 
DO Collaboration performed a search for pair production 
of leptoquarks decaying into a quark and a neutrino 14( , 
which leads to the same topology. The analysis reported 
here is largely based on that leptoquark search. In the 
following, only a brief summary of the common aspects 
is given, while the specific features relevant for the stop 
search are presented in greater detail. The main differ- 
ences arise from the LSP mass, which leads to smaller jet 
transverse energies and to a reduced $t, compared to the 
case of leptoquark decays which involve nearly massless 
neutrinos. Another characteristic feature of stop decays 
is that charm jets are produced, while first-generation 
leptoquarks decay to light-flavor jets. 



A thorough description of the DO detector can be found 
in Ref. [lij]. The central tracking system consists of a 
silicon microstrip tracker and a fiber tracker, both lo- 
cated within a 2T superconducting solenoidal magnet. 
A liquid-argon and uranium calorimeter covers pseudo- 
rapidities \rj\ < 4.2, where rj = — In [tan (0/2)] and 8 is 
the polar angle with respect to the proton beam direc- 
tion. An outer muon system, covering \rj\ < 2, consists 
of layers of tracking detectors and scintillation counters 
on both sides of 1.8 T iron toroids. 

For this search, 14 million events collected from 
April 2003 to August 2004 with a jets + trigger were 
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analyzed, corresponding to an integrated luminosity 1 of 
360 pb -1 . The offline analysis utilized jets reconstructed 
with the iterative midpoint cone algorithm 17| with a 
cone size of 0.5. Only jets with transverse momentum 
Pt > f 5 GeV were considered in the analysis. The fir 
was calculated using all calorimeter cells, corrected for 
the energy calibration of reconstructed jets, as deter- 
mined from the transverse momentum balance in pho- 
ton+jet events, and for the momentum of reconstructed 
muons. 

Signal efficiencies and SM backgrounds were evaluated 
using a full GEANT-3 18| based simulation of events, with 
a Poisson average of 0.8 minimum-bias events superim- 
posed, corresponding to the luminosity profile of the data 
sample analyzed. These simulated events were recon- 
structed in the same way as the data. In the bulk of 
events from QCD multijet production, no significant fir 
is expected. Jet energy mismeasurements due to the lim- 
ited detector resolution may however lead to large mea- 
sured fir values. This "instrumental background" was 
not simulated, and its contribution estimated directly 
from the data. In the following, "standard model (SM) 
background" stands for "non-QCD standard model (SM) 
background." Leptonic W decays, as well as Z — > vv are 
sources of energetic neutrinos, hence of genuine fir- The 
SM processes expected to yield the largest background 
contributions are therefore vector boson production in 
association with jets. They were generated with ALP- 
GEN 1.30, interfaced with pythia 6.202 for the 
simulation of initial and final state radiation and for jet 
hadronization. The PDFs used were CTEQ5L[21]. The 
NLO cross sections for vector boson production in asso- 
ciation with jets were calculated with mcfm 3.4.4 
Vector-boson pair, tt, and single top quark production 
were also considered. Signal samples of 10 000 events 
were generated with pythia and the CTEQ5L PDFs 
for stop masses ranging from 95 to 145 GeV and for Xi 
masses from 40 to 70 GeV, both in steps of 5 GeV. 

The following selection criteria were applied, indepen- 
dent of the stop and Xi masses: there had to be at least 
two jets; the vector sum fir of all jet transverse momenta 
(I/It = | Sjcts'PVI) as wen as the missing transverse en- 
ergy had to exceed 40 GeV; the leading and subleading 
jets (where jets are ordered according to their transverse 
momentum) had to be central (|f?det| < 1-5, where 77det 
is the pseudorapidity measured from the detector cen- 
ter), with transverse momenta exceeding 40 and 20 GeV, 
respectively, and they had to be confirmed by charged 
particle tracks [3] ; the acoplanarity A<I> of the two lead- 
ing jets had to be smaller than 165°, where A$ is the 
difference between the two jet azimuthal angles; the lon- 



This value differs from the one used in Ref. |14| due to a recent 
adjustment of the DO luminosity constant Hal . 



gitudinal position of the primary vertex had to be less 
than 60 cm away from the center of the detector. At this 
point, 99 884 events were selected, largely dominated by 
instrumental background from multijet events. The ef- 
ficiency for a reference signal with mi = 140 GeV and 
m^o = 60 GeV was 30%. 

The jet multiplicity distribution revealed that most 
of the selected events contained at least three jets, due 
to the acoplanarity requirement. Therefore, only events 
containing exactly two jets were retained, leaving 27 853 
data events with an efficiency of 22% for the reference 
signal. The inefficiency associated with the rejection of 
events with more than two jets was evaluated, based on 
studies of jet multiplicities in real and simulated Z — > ee 
events with at least two jets, where the two leading jets 
fulfilled similar selection criteria as in the analysis. This 
study also showed that the kinematic variables used in 
the analysis were adequately simulated. Standard model 
backgrounds from W — > iu+jet processes were greatly 
reduced by requiring that there be no isolated electron 
or muon with pr > 10 GeV, and no isolated charged par- 
ticle track with p T > 5 GeV [3 . This retained 22106 
data events, with an efficiency of 19% for the reference 
signal. 

Most of the remaining instrumental background was 
eliminated by the following requirements. The fir had 
to exceed 60 GeV, and the difference V — A$ max — A$ min 
had to be smaller than 120°, where A$ m i n and A$ max 
are the minimum and maximum of the azimuthal angles 
between the fir direction and the directions of the two 
jets, respectively. These criteria take advantage of the 
facts that, for the instrumental background, the fir dis- 
tribution is steeply decreasing, and its direction tends 
to be close to that of a mismeasured jet. In addition, 
the asymmetry A — (fir ~ fir) /(fir + fir) was re- 
quired to be larger than —0.05. This variable is sensi- 
tive to the amount of energy deposited in the calorime- 
ter that was not clustered into jets. It can be seen in 
Fig.[T] that both T> and A are effective in discriminating 
SM backgrounds and signal from the instrumental back- 
ground. After these requirements, 1 348 data events were 
retained, while 1 292 ± 45 events were expected from SM 
backgrounds, where the uncertainty is statistical. The 
efficiency for the reference signal was 13%. There was 
no evidence at this point for any significant instrumental 
background remaining. This background has therefore 
been neglected in the following. 

To increase the search sensitivity, advantage was then 
taken of the presence of charm jets in the signal. A 
lifetime-based heavy-flavor tagging algorithm was used 
for this purpose, which involves a probability built from 
the impact parameter significances of the tracks belong- 
ing to a jet [23|]. The impact parameter of a track is 
its distance of closest approach to the event vertex, in 
a plane perpendicular to the beam axis, and the signifi- 
cance is obtained by normalization to the impact param- 
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FIG. 1: Distributions of the asymmetry A = ($t — ]^t)/(^t + $t) with the cut on V = A$ max — A$ m i n inverted (top-left) 
or applied (bottom-left) and of T> with the cut on A inverted (top-right) or applied (bottom-right) for data (points with error 
bars), for SM backgrounds (filled histogram), and for a signal with = 140 GeV and m^o = 60 GeV (hatched histogram). The 
JpT cut at 60 GeV has been applied. In the bottom plots, the excesses in data for A < —0.05 and for T> > 120° are attributed 
to the residual non-simulated instrumental background. 



eter uncertainty. This probability is constructed such 
that its distribution is uniform for light-flavor jets and 
peaks towards zero for heavy-flavor jets. In order to 
cope with differences in track reconstruction efficiencies 
in data and in simulation, the heavy-flavor tagging algo- 
rithm was applied directly only to the data, while flavor- 
dependent tagging probabilities measured in dedicated 
data samples were applied to the simulated jets. The 
probability cut used in this analysis was such that typi- 
cally 4% of the light-flavor jets were tagged (central jets 
with pt s* 50 GeV). The corresponding typical tagging 
efficiencies for c and b quark jets were 30% and 65%, re- 
spectively. Jets resulting from r decays were tagged with 
a typical efficiency of 20%. By requiring that at least 
one jet be tagged, 183 data events were selected, while 



186±16 SM background events were expected, where the 
uncertainty is statistical. The efficiency for the reference 
signal was 6.5%. 

Since the signal topology depends on the stop and x? 
masses, additional selection criteria on three kinematic 
variables were simultaneously optimized for each mass 
combination. These variables were the scalar sum Hj 1 — 
Sjets Y^t\ °f the jet transverse momenta in steps of 
20 GeV, i? T in steps of 10 GeV, and S = A$ max + A$ min 
in steps of 10°. It can be seen in Fig. [2] that this last 
variable provides good discrimination between signal and 
SM backgrounds. For Ht and Ifix, the selection retained 
events above the cut value, while for S, events below the 
cut value were selected. For each stop and \i mass com- 
bination tested, all sets of cuts were considered. For each 
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set, the value (CL S ) of the signal confidence level 24| ex- 
pected if only background were present was computed, 
with the systematic uncertainties discussed below taken 
into account. For a given stop mass, the expected lower 
limit on rra^o was determined as the Xi mass for which 
(CL S ) — 5%, by interpolation across the m^o values 
tested. The set leading to the largest expected lower limit 
on m^o was selected as the optimal one for the stop mass 
considered. In all IpT cut at 60 GeV was selected. 

The results of the optimization for the other variables are 
given in TableHJ together with the numbers of events se- 
lected in the data and expected from SM backgrounds. 
Signal efficiencies and numbers of signal events expected 
are given m TableHJfor three mass combinations close to 
the edge of the sensitivity domain of the analysis. 

The distribution of Ht shown in Fig. [5] and the final 
distribution of IpT shown in Fig.[3]were obtained after op- 
timization for a stop mass of 140 GeV. An excess at large 
JpT is observed in the data with respect to the expec- 
tation: there are eight data events with Ifr > 150 GeV, 
while 3.2±1.4 events are expected from SM backgrounds. 
A detailed scrutiny of those events was performed, that 
did not reveal any anomaly such as clusterings in some 
of the kinematic variables, signs of leptons unidentified 
by the standard algorithms, heavy flavor tagging proba- 
bilities different from what is observed in the rest of the 
selected events. The data taking conditions were also 
carefully checked for signs of detector malfunctions and 
visual scans were performed. It can also be noted that 
such large I$t values are beyond what is expected from 
a stop signal. 



TABLE I: Results of the optimization: stop mass range in 
GeV, Ht cut value in GeV, and S cut value in degrees. In 
all cases, a IpT cut at 60 GeV was selected. The numbers of 
events observed and expected from SM backgrounds are also 
given; the first uncertainties are statistical, and the second 
systematic. 





Ht 


S 


# observed 


# expected 


95 - 115 


> 80 


< 260 


68 


59.9 ± 9.6 tl 1 ? 


120 


> 80 


< 280 


89 


86.4 ± 11.3 ±il;| 


125 - 140 


> 120 


< 280 


50 


47.0 ±8.0 tjl 


145 


> 120 


< 300 


57 


53.8 ± 8.3 l9° 2 8 



The SM background composition is detailed in Ta- 
blc lIIII for the selection optimized for = 140 GeV. As 
expected, the largest contributions come from (Z — > vv 
and W — > ^)±light-fiavor jets. This is due to the loose 
heavy-flavor tagging criterion which was selected in order 
to be efficient for charm jets. Vector boson production 
with heavy-flavor jets gives rather small contributions be- 
cause of the comparatively small cross sections. 

Systematic uncertainties were evaluated for each com- 
bination of stop and Xi masses, according to the cor- 
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FIG. 2: Distributions of S = A$ max + A$ m i n before op- 
timization (top), and of Ht after optimization for m t - = 
140 GeV but with the cut on Ht removed (bottom) , for 
data (points with error bars), for SM backgrounds (filled his- 
togram) , and for a signal with m t - = 140 GeV and m^o = 
60 GeV (hatched histogram). 



TABLE II: For three stop and x? mass combinations, in GeV, 
signal efficiencies (Eff.) and numbers of signal events ex- 
pected, where the first uncertainties are statistical and the 
second systematic. The stop pair production cross section 
upper limits at 95% C.L. are also given (<7ul), as well as the 
NLO theoretical cross section (otii), both in pb. 



(m t -,m s o) 


Eff. (%) 


# expected 






(100,55) 


0.75 


40.4 ±4.6 tH 


15.8 


15.0 


(120,65) 


2.04 


40.0 ±2.8 


6.57 


5.43 


(140,60) 


3.74 


30.3 ±1.6 tH 


2.38 


2.25 
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FIG. 3: Final 1$t distribution for data (points with error 
bars), for SM backgrounds (filled histogram), and, on top of 
the SM backgrounds, for a signal with nif — 140 GeV and 
m^o = 60 GeV (hatched histogram). 



TABLE III: Numbers of events expected from the various SM 
background processes in the selection optimized for = 140 
GeV. The uncertainties are statistical. In the vector boson ± 
jets backgrounds, "jet" stands for "light-flavor jet." 



SM process 


# expected 


Z —* vv+]ets 


13.9 ± 2.8 


Z — > vv+cc 


1.7 ±0.4 


Z -> vv+bb 


3.5 ±0.2 


W -> &/+jets 


19.5 ±7.4 


W — > £u+(cc or c+jet) 


1.8 ±0.5 


w -> eis+bh 


1.5 ±0.2 


tt and single top 


4.1 ±0.2 


WW, wz, zz 


1.1 ±0.2 


Total 


47.0 ± 8.0 



responding optimized selection criteria. They are listed 
below for the reference signal. The following are fully cor- 
related between SM-background and signal expectations: 
from the jet energy calibration and resolution, L"g 3 % for 
the SM background and ^4% for the signal; from the jet 
multiplicity cut, 3%; from the trigger efficiency, 2% after 
all selection cuts; from the heavy-flavor tagging, 6% for 
the SM background and 7% for the signal; from the inte- 
grated luminosity of the analysis sample, 6%. In addition 
to the 17% statistical uncertainty of the simulation, the 
normalization of the SM background expectation carries 
a 13% uncertainty, as inferred from a comparison of data 
and simulated (Z — > ee) + 2-jet events. The statistical 
uncertainty of the signal simulation is 5%. Finally, the 
uncertainty on the signal efficiency due to the PDF choice 
was determined to be t.4%, using the CTEQ6.1M error 
set [HI. 



As can be seen in TableU no significant excess of data 
was observed in any of the optimized selections. Signal 
production cross section upper limits were therefore de- 
rived with the above systematic uncertainties taken into 
account. Examples are given in TablellTl together with 
the corresponding theoretical cross sections. To deter- 
mine an exclusion domain in the (m^m^o) plane, the 
following procedure was used. For a given the signal 
confidence level CL S was computed as a function of m^o 
in the modified frequentist approach 24], and the 95% 
C.L. lower limit on m^o was determined as the Xi mass 
for which CL S = 5%. In this procedure, the theoreti- 
cal NLO cross sections predicted by prospino with the 
CTEQ6.1M PDFs were used. The nominal cross section 
was obtained for ji r t — mj. Theoretical uncertainties 
on the stop pair production cross section arise from the 
choices of PDFs and of renormalization and factorization 
scale. The variations observed with the CTEQ6.1M er- 
ror PDF set, as well as the changes induced when /j, r f is 
modified by a factor of two up or down, result in a typi- 
cally ±20% change in the theoretical cross section when 
combined in quadrature. The exclusion contour in the 
(mj,m^o) plane thus obtained is shown as a solid curve 
in Fig. [4] for the nominal production cross section. The 
corresponding expected exclusion contour is shown as a 
dashed curve. The effect of the PDF and scale uncer- 
tainties on the observed exclusion contour is shown as a 
shaded band. 



This analysis, performed under the assumption that 
the stop decays exclusively into a charm quark and the 
lightest neutralino, extends the stop and y? mass domain 
excluded by previous experiments [ll|, [lj, [HI . For the 
nominal stop pair production cross section, the largest 
stop mass excluded is 141 GeV, obtained for to^o = — 
irib ~~ m 'W = 55 GeV. Taking into account the theoretical 
uncertainty on the production cross section, the largest 
stop mass limit is 134 GeV, obtained for m^o = 48 GeV. 
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FIG. 4: Domain in the (m^m^o) plane excluded at the 95% 
C.L. by the present search (region below the solid curve), un- 
der the assumption that the stop decays exclusively into cx? 
and for the nominal production cross section. The expected 
exclusion contour is shown as a dashed curve. The effect of 
increasing or decreasing the production cross section by its 
uncertainty due to the PDF and fj, r f choices is indicated for 
the observed exclusion contour by the shaded band. Results 
from previous searches for stop pair pro duction in the i — > ex? 
decay channel are also indicated [lll.[l2l . The dark shaded 
band at small m t - — m^o is excluded by Ref. [25|. The LEP 
results are shown for two values of 8, the mixing angle in the 
stop sector. 



[f] Visitor from Helsinki Institute of Physics, Helsinki, Fin- 
land. 

[1] H.E. Haber and G.L. Kane, Phys. Rep. 117, 75 (1985). 
[2] See for instance V. Barger, M.S. Berger and P. Ohmann, 
Phys. Rev. D 49, 4908 (1994). 



[3] 
[4] 

[5. 
[6 

[<: 

[8 

[9 
[10 

[11 



[12 

[13 
[14 

[is; 
[ie; 

[17; 

[is; 

[19 
[20 

[21 

[22 

[23 
[24 

[25; 



J. Ellis and S. Rudaz, Phys. Lett. B 128, 248 (1983). 
M. Quiros, Nucl. Phys. Proc. Suppl. 101, 401 (2001), and 
references therein. 

P. Fayet, Phys. Lett. B 69, 489 (1977). 

J. Ellis et al, Nucl. Phys. B 238, 453 (1984). 

K.I. Hikasa and M. Kobayashi, Phys. Rev. D 36, 724 

(1987). 

C. Boehm, A. Djouadi and Y. Mambrini, Phys. Rev. D 
61, 095006 (2000). 

W. Beenakker et al, Nucl. Phys. B 515, 3 (1998). 

J. Pumplin et al., J. High Energy Phys. 0207, 012 (2002); 

D. Stump et al, ibid., 0310, 046 (2003). 
LEPSUSYWG, ALEPH, DELPHI, L3 and OPAL Col- 
laborations, note LEPSUSYWG/04-02.1 



| http://lepsusy.web.cern.ch/lepsusy/Welcome.html |. 
T. Affolder et al. (CDF Collaboration), Phys. Rev. Lett. 
84, 5704 (2000). 

V.M. Abazov et al. (DO Collaboration), Phys. Rev. Lett. 
93, 011801 (2004). 

V.M. Abazov et al. (DO Collaboration), Phys. Lett. B 
640, 230 (2006). 

V.M. Abazov et al. (DO Collaboration), Nucl. Instrum. 
and Methods A 565, 463 (2006). 

T. Andeen et al, "Adjustments to the DO Experi- 
ment's Measured Luminosity for Tevatron Run Ha," 
FERMILAB-TM-2365, in preparation. 

G. C. Blazey et al, "Run II Jet Physics: Proceedings of 
the Run II QCD and WeaJc Boson Physics Workshop," 
arXiv: |hep-ex/ 0005012 

R. Brun and F. Carminati, CERN Program Library Long 

Writeup W5013, 1993 (unpublished). 

M.L. Mangano et al, J. High Energy Phys. 0307, 001 

(2003). 

T. Sjostrand et al, Comput. Phys. Commun. 135, 238 
(2001). 

H. L. Lai et al. Eur. Phys. J. C 12, 375 (2000). 

J. Campbell and R.K. Ellis, Phys. Rev. D 60, 113006 
(1999). 

B. Clement, Ph.D. thesis, Universite Louis Pasteur, 

Strasbourg, N° d'ordre IPHC 06-004, N° d'ordre ULP 

5086, FERMILAB-THESIS-2006-06 (2006). 

T. Junk, Nucl. Instrum. and Methods A 434, 435 (1999); 

A. Read, in "First Workshop on Confidence Limits," 

CERN Report No. CERN-2000-005, 2000. 

A. Heister et al. (ALEPH Collaboration), Phys. Lett. B 

537, 5 (2002). 



